Observable Vortex Properties in Finite Temperature Bose Gases 



(N 

O 

(N 

> 

o 

O: 



C/2 



o 

> 
m 

d 



X: 



A. J. AUcn^, E. Zarcmba'^, C. F. Barcnghi^, and N. P. Proukakis^ 
^ Joint Quantum Centre (JQC) Durham-Newcastle, School of Mathematics and Statistics, 
Newcastle University, Newcastle upon Tyne NEl 7RU, England, UK. 
^Department of Physics, Engineering Physics and Astronomy, 
Queen's University, Kingston, Ontario K7L 3N6, Canada. 
(Dated: November 1, 2012) 

We study the dynamics of vortices in finite temperature atomic Bose-Einstein condensates, 
focussing on decay rates, precession frequencies and core briglitness, motivated by a recent 
experiment (Freilich et al. Science 329, 1182 (2010)) in which real-time dynamics of a single vortex 
was observed. Using the ZNG formalism based on a dissipative Gross-Pitaevskii equation for the 
condensate coupled to a semi-classical Boltzmann equation for the thermal cloud, we find a rapid 
nonlinear increase of both the decay rate and precession frequency with increasing temperatures. 
The increase, which is dominated by the dynamical condensate-thermal coupling is also dependent 
on the intrinsic thermal cloud coUisional dynamics; the precession frequency also varies with the 
initial radial coordinate. The integrated thermal cloud density in the vortex core is for the most 
part independent of the position of the vortex (except when it is near the condensate edge) with its 
value increasing with temperature. This could potentially be used as a variant to the method 
of Coddington et al. (Phys. Rev. A 70, 063607 (2004)) for experimentally determining the 
temperature. 

PACS numbers: 03.75.Kk, 67.85.-d 



I. INTRODUCTION 



The existence of quantised vortices in Bosc-Einstein 
condensates (BECs) is a feature of the superfluid na- 
ture of these systems [l|. Since the first experimental 
realisation of vortices in BECs 043 j there have been 
numerous investigations into the nature of their dynamics 
(see e.g. la), ranging from their nucleation as a result of 
rotation [6[ (see also and references therein), to the 
effect of the inhomogeneous density on the velocity of 
a vortex to the decay of a single vortex induced 

by sound emission and the formation and struc- 

ture of vortices in multicomponent [l6l - [l9| and spinor 
BECs [20l - [2^ . Vortices are now routinely created in a 
variety of waro, including stirring the condensate using a 
laser beam |23| . j letting a soliton decay via the snake 
instability j25| . phase imprinting [26j . and by a rapid 
quench through the transition temperature for the onset 
of Bose-Einstein condensation (i.e. the Kibblc-Zurek 
mechanism) [13, . These techniques have opened the 
door to the study of more complicated configurations, 
for instance, the motion of vortex dipoles l28l. [29| . the 
formation of multiply charged vortices [26l. ISOf. vortex 
lattices H, 0, [Sll, |32| and, more recently, the creation 
of a small tangle of vortices [33l - [36| . The lifetime of 
vortex structures created in the laboratory can range up 
to several seconds. The reason for their eventual demise 
is thought to be thermal dissipation [s^ . [37| . At finite 
temperature, atom-atom interactions cause the vortex to 
lose energy; this results in the vortex spiralling out of a 
harmonically trapped condensate (see Fig.[3J top.); thus, 
the lifetime of a vortex is severely reduced with increasing 
temperature (38| . Earlier work on finite temperature 
vortex dynamics has confirmed this effect p38l - l46| . Along 



with the decay, the precession frequency of a singly- 
charged, harmonically trapped vortex has also been 
found to vary with temperature [H, IH, [i^l ■ 

The motivation behind our work is twofold. Firstly we 
revisit the problem of vortex decay rate and precession 
frequency, using the only implementable model to date 
which includes the full thermal cloud dynamics [48| : 
other approaches typically only include the dynamics 
up to a cutoff within classical field theory (which is 
more suitable in the fluctuation-dominated regime, very 
close to Tc - see, for example, [49|V Our study is 
motivated by a recent experiment [28| which followed 
the motion of a single vortex in a harmonically trapped 
condensate using a new imaging technique. Secondly, 
we are interested in revisiting the experimental proposal 
of Coddington et al. [s^l regarding the use of finite 
temperature effects on vortex core 'brightness' as a 
possible tool for thermometry; at low temperatures it 
is difficult to accurately extract the temperature of a 
system of bosons, and since the mean-field felt by the 
thermal cloud in the region of the vortex core is much 
less than anywhere else in the condensate, the vortex core 
was said to act as a 'thermal-atom concentration pit' [s^] ■ 
Virtanen et al. [5l| have theoretically investigated the 
core brightness over a small range of temperatures, for a 
central vortex using the Hartree-Fock Bogoliubov theory 
within the 'Popov' approximation (52j . and found it to 
increase with temperature, as expected. We revisit this 
problem, using a fully dynamical theory (isl . [53j . and 
determine the dependence of core contrast on both the 
radial component of the vortex and the temperature. 

Our approach is based on the formalism of Zaremba, 
Nikuni and Griffin (ZNG) 0, HI, where the Gross- 
Pitaevskii equation (GPE) is generalised by the inclusion 
of the thermal cloud mean field, and a dissipative/source 



2 



term which is associated with a colhsion term in a 
semi-classical Boltzmann equation for the thermal cloud. 
This approach has previously been used to successfully 
describe the damping of condensate collective modes [5J- 
[ssj and macroscopic excitations [H, HI, [g^ in the mean- 
field dominated regime at finite temperature; the method 
reduces to the damped two-fluid equations of ^Hc in the 
hydrodynamic limit [48. i51. JiSj . This is the only model 
implemented to date which self-consistently accounts for 
all collisional dynamics of the system and is suited to 
elevated temperatures excluding the region of critical 
fluctuations. 

The plan of this article is the following. In Section [IT] 
we briefly review the ZNG equations (with some further 
details in Appendix |^ and apply them to the problem 
of vortex decay at finite temperatures. Following the 
procedure of Jackson et al. [38|, we extract a decay 
rate for various temperatures for a fixed number of 
condensate atoms, and analyse the effect of collisions 
on the decay rate. In Section IIIIl we use the same 
parameters to assess how increasing temperature affects 
the precession frequency of a vortex. In particular in 
Sec nil Al we extract precession frequencies of vortices 
for the parameters of the Freilich et al. [2^ experiment, 
where the total number of atoms are fixed. In Section ITVl 
we investigate vortex core brightness, and how it changes 
with radial position and temperature, and conclude by 
summarising our findings in Section FVl 



II. BACKGROUND, THEORY AND 
MOTIVATION 



In typical vortex experiments, the vortex core size is 
smaller than the wavelength of the light used to image it. 
As a result, to make the vortex visible it is necessary to 
expand the condensate [2^. However, as a consequence 
of this expansion, the condensate is destroyed, requiring 
successive reproducible runs to observe time evolution of 
the vortex. 

Freilich et al. [l^ developed an imaging technique, 
which involves the repeated extraction and expansion of 
approximately 5% of the condensate atoms, thus enabling 
a scries of images of the same condensate containing 
the vortex to be created. This technique allowed the 
precession frequency of the vortex to be measured. 

We will now briefly discuss the ZNG equations and the 
advantages of using such a model before presenting our 
results. 



A. The ZNG formalism 

An extensive review of the ZNG formalism and its 
derivation can be found elsewhere [H,!!!] . Here it suffices 
to briefly outline the methodology. The formalism is 
based on the following closed set of equations (where 
the explicit dependence on r and t is suppressed for 



convenience 
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i-n— = — - — 
dt \ 2m 



+ T4xt + K + 2?i] - ii? U , (1) 



a/ 

dt 



+ - • Vr/ - (Vrt/cff ) • (Vp/) = Ci2 [/, 0] + C22 [/] . 



(2) 



Here, (j) is the condensate wavefunction and / is a 
phase-space distribution function. Eq. ([T]) is a finite 
temperature generalisation of the usual Gross-Pitaevskii 
equation which is modified by the addition of a thermal 
cloud mean- field potential 2gn, and a dissipativc/source 
term — i_R. It is coupled to the quantum Boltzmann 
equation (QBE, Eq. ©) for the thermal cloud phase- 
space distribution. The condensate density is written 
as Tic — Kxt(r) is the external trapping poten- 

tial, and the interaction strength between the atoms 
is g — ATrh^as/m, where as is the s-wave scattering 
length and m is the atomic mass. The thermal cloud 
density is recovered from the phase-space distribution 
function via an integration over all momenta, n(r, t) = 
/ dp/(27r?i)3/(p,r,t) and U^s = K=xt(r) + 2.gK(r,t) + 
?i(r,^)] is the mean-field potential acting on the thermal 
atoms. The quantities C22 and C12 are collision integrals 
(definitions of these terms can be found in Appendix . 
C22 describes the redistribution of thermal atoms as 
a result of collisions between two thermal atoms, i.e. 
the usual Boltzmann equation collision integral, while 
C12, which is closely related to the dissipation term 
iR, describes the change in the phase-space distribu- 
tion function /(p, r, t) as a result of particle-exchanging 
thermal atom-condensate collisions. Although other 
methods have been put forward for finite temperature 
vortex dynamics, to the best of our knowledge, this is 
the only method that self-consistently accounts for the 
redistribution of thermal particles and its effect this has 
on the condensate dynamics [69| . 



B. Vortex Decay 

Equations dH) and ^ have previously been used to 
study finite temperature vortex dynamics by Jackson 
et al. d^. In that work, the authors show that the 
decay rate of a vortex increases rapidly with increasing 
temperature. In these simulations, the system size was 
fixed to a constant total number of atoms, N^qt = 
10, 000; by increasing the temperature of the system, the 
number of thermal atoms increased and consequently the 
condensate size decreased. As both of these variations 
affect the vortex dynamics, in order to isolate the effect of 
increasing thermal atom number, in this work we instead 
initially perform simulations at different temperatures 
for a fixed number of condensate atoms, Nc — 10, 000 
^^Rb atoms. By fixing the number of condensate 
atoms, the total number of atoms, iVTOT increases with 



3 



temperature, T. As a result, the critical temperature, T^, 
for the onset of Bose-Einstein condensation is a function 
of T in these simulations. We estimate T'c(-/Vtot) by 
means of the ideal gas expression (note that, for A^tot = 
10, 000 atoms, the T = value of the critical temperature 
is Tc = 177nK). 

The geometry is a fully three-dimensional (3D) har- 
monic trap, Kxt(r) = m/2(a;^(x^ 4- y"^) + uj^z^) with 
trapping freq uencies, uj± ^ 2tt x 129Hz and = 
■v/Swj^Hz [3g. The purpose of the significantly tighter 
trapping frequency in the axial direction is to ensure that 
the vortex remains relatively straight along its length, as 
shown in Fig. [T] (left). This figure contains 3D isosurface 
plots of the condensate (left) and thermal cloud (right) 
densities. We have chosen a high density surface of the 
thermal cloud in order to show how the noncondensate 
atoms fill the vortex core and concentrate around the 
edges of the condensate. Figure [2l shows the densities 




FIG. 1. (Color online) 3D isosurface plots of low condensate 
density (left, red) and high thermal cloud density (right, blue) 
for a cloud of Nc — 10, 000 ^'^Rb atoms at a temperature 
of 0.5Tc for the trapping parameters as described in the 
text. Notice the tubular isosurface of the thermal cloud at a 
position corresponding to the vortex core in the condensate. 

of both components integrated along the z— direction. 
The relatively flat profile of integrated thermal density 
in Fig. [2] (right), illustrates how the thermal cloud 
surrounds the condensate, with the regions of highest 
concentration of noncondensate atoms corresponding to 
the condensate edge as well as position of the vortex 
(these are the areas where the thermal cloud feels the 
lowest mean- field repulsion from the condensate). 

We extract the decay rate, 7, for vortices of different 
initial positions, by fitting the vortex radial variable to an 
exponential, ry{t) = rge'''*, where vq is the initial vortex 




FIG. 2. (Color online) Two-dimensional (2D) densities, 
integrated along the direction for the condensate (left) and 
thermal cloud (right). Notice the effect of the vortex core 
on the thermal cloud. Trapping parameters as in text for: 
T/Tc = 0.5, iVc = 10, 000 **^Rb atoms. 



position. Results corresponding to three temperatures 
are plotted in Fig. [31 where {xy,yy) trajectories are also 
shown. We express the vortex radial coordinate in terms 
of the Thomas-Fermi radius which is defined as i?TF = 
y^2fj,/m(jjj_, where n is the chemical potential. 
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FIG. 3. (Color online) Top: {xv,yv) trajectory of vortex 
initially at ro — 0.4i?TF at temperatures T/Tc — 0.3 (53nK) 
(left, black), 0.5 (89nK) (middle, red) and 0.6 (124nK) (right, 
blue). These trajectories have been smoothed to remove 
any numerical 'jitter' arising as a result of our tracking 
routine. Bottom: Corresponding values of decay rate, 7, 
for vortices of variable initial position, ro = + y^, at 

these temperatures. The filled symbols indicate the position 
of the vortex subsequently analysed in Fig. [l] Trapping 
parameters as in text for Nc — 10, 000 *^Rb atoms, Tc varies 
between (175 — 205)nK for the range of temperatures shown. 
iiTF — 5l± with the maximum values of the thermal cloud 
density occurring at ~ 0.85i?TF in all cases. 

Increasing the temperature from 0.3Tc to 0.7rc in- 
creases the decay rate significantly, in agreement with 
previous findings [H, [H, Hi], [t^]- It is also apparent 
that at all temperatures, moving the vortex initial radial 
coordinate rp away from the centre and closer to the 
condensate edge {Rtf — 5/_l) results in an increase 
in vortex decay rate, 7. This feature is notably more 
pronounced at the higher temperatures, and is caused by 
the higher thermal cloud density towards the edge of the 
condensate (the maximum thermal cloud density occurs 
at the position w 0.85i?TF) which arises because of the 
weaker mean-field repulsion. We now discuss the effect of 
the individual collision terms of the QBE, Eq. ©, on the 
vortex decay rate as well as on the precession frequency 
of a vortex. We measure the precession frequency of a 
vortex as the inverse time it takes for one oscillation of 
the trap, calculated by averaging the oscillation time over 
the first three oscillations. 
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1. Effect of collisions on the vortex decay rate 

We can carry out numerical simulations with the full 
QBE of Eq. however, we can also simulate the effect 
of the thermal cloud with a combination of these collision 
terms, or without their inclusion at all. 

In this way, we can determine the contributions to 
the vortex decay rate coming from different collisional 
processes. Fig. |4] gives the results of vortex decay rate, 
7, and precession frequency, uj^ for a vortex having an 
initial radial offset from the trap centre of tq ~ 0.26i?TF, 
for various temperatures. The filled circles show the 
results for the full QBE simulation. Results are also 
shown for when the thermal-thermal, C22, collisions are 
neglected but the particle-exchanging, C12, collisions are 
included (magenta squares), and vice versa (blue pluses, 
consequently, when C12 = the dissipative/source term 
iR will also be zero). In the case when both of these 
terms are neglected (open, green circles), the QBE is 
propagated in time according to free streaming terms 
(left hand side of Eq. © ) , which ensures that the value of 
the phase-space distribution function remains the same 
along a trajectory in phase-space. We have also obtained 
results for a static thermal cloud (red stars), in which 
the mean-field of the equilibrium thermal cloud density 
(obtained in the absence of a vortex) is included in the 
solution of the GPE. These results highlight the crucial 
role of all collision processes in determining the actual 
decay rate and precession frequency of a vortex. 

Focussing initially on the decay rate (Fig. |31 top), 
we see that the largest contribution to vortex decay 
comes from the particle exchanging, C12, collisions. 
When the C12 collisions are not included, the decay 
rate is reduced significantly. The thermal-thermal, C22, 
collisions have a noticeable effect only when they arc 
included together with the particle-exchanging, C12, 
collisions. This implies that C22 collisions affect the 
decay rate indirectly through a modification of the C12 
collision rates. 

For the precession frequency of a vortex (Fig. IH 
bottom), we see increasing values with increasing tem- 
perature, with the effect of the inclusion/absence of 
the different collision terms being more apparent at 
the higher temperatures. The largest influence on 
the precession frequency again arises from the particle- 
exchanging, C12, collisions - an effect which is intuitive 
since these collisions cause the vortex to lose energy and 
move out of the condensate radially. 

To summarise, our analysis demonstrates the impor- 
tance of including the full dynamics of the thermal 
cloud, i.e. all of the collision terms of the QBE when 
modelling vortex dynamics. All further results quoted 
have been simulated with all collision terms included 
in the propagation of Eq. In the next section, we 
analyse further the dependence of the vortex precession 
frequency on temperature, as well as with increasing 
vortex radial coordinate. 
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FIG. 4. (Color online) Vortex decay dynamics as a function 
of temperature: values of vortex decay rate, 7, (top) and 
corresponding vortex precession frequencies, ujv, (bottom). 
Results for different levels of approximation are indicated by 
(i) solid, black dots: all collision (Ci2,C22) processes, (ii) 
open, magenta boxes: particle-transferring (C12) collisions 
only, (iii) blue crosses: thermal-thermal (C22) collisions only, 
(iv) open, green circles: no collisions, and (v) red stars: 
static thermal cloud approximation. An initial radial vortex 
offset of ro ~ 0.26-Rtf, as highlighted in Fig. [3] is used for 
all simulations. Clearly all collision mechanisms contribute 
significantly to the decay. Trapping parameters as in text 
for Nc = 10,000 ^■'Rb atoms with varying between 
(175 — 280)nK for the results shown here. 



III. PRECESSION FREQUENCY OF A VORTEX 

Our motivation for investigating further the effect of 
finite temperature on vortex precession arises from the 
development of a novel imaging technique [28j . enabling 
real-time vortex dynamics to be observed. To first 
understand the effect of the initial vortex position on 
its subsequent precession, we begin our analysis by 
extracting the precession frequencies for vortices having 
various different initial positions for a fixed temperature 
(Fig.[5]inset). The simulations are again carried out using 
the trapping parameters quoted in Sec. IIIBI and a fixed 
number of condensate atoms, Nc = 10, 000 for various 
temperatures. 

In the inset of Fig. [5] we plot the precession fre- 
quency as a function of vortex radial coordinate, r^{t) = 
■\/a;t,(t)^ -f yy{t)'^, for initial positions rp in the range 
(0.1 - 0.8)i?TF at the temperature T/Tc = 0.6. 

The general increase of the precession frequency with 
increasing vortex radial coordinate, r^, is apparent. The 
fact that all the points appear to lie on a common 
curve confirms numerically an important anticipated 
feature of vortex dynamics: the frequency of vortex 
precession depends on the instantaneous radial position 
and not on the history of how the vortex arrived at that 
position. We can therefore average over all the points 
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FIG. 5. (Color online) Inset: Frequencies of vortices with 
different initial positions ro in the range (0.1 — 0.8)-Rtp 
(shown by dots of different colours) at a temperature of 
T/Tc = 0.6. Main: The corresponding averaged curve, for this 
temperature, T/Tc = 0.6 (124nK) (blue dots and connecting 
line). For the remaining temperatures, from bottom to top 
with increasing temperature, T/Tc = (green diamonds), 
0.3 (53nK) (black triangles), and 0.5 (88.5nK) (red squares), 
precession frequencies have been extracted as a function of 
radial coordinate. Trapping parameters as in text for Nc ~ 
10, 000 *'^Rb atoms, -Rtf ~ 5l± for these temperatures and 
Tc varies between (175 — 205)nK for the results shown here. 



to generate a representative curve for the frequency as 
a function of radial coordinate. This curve is shown in 
the main part of Fig. [5] (bkie dots). Similar simulations 
were performed for the other temperatures displayed in 
Fig. O Fewer points are shown as these fully capture the 
investigated behaviour of the frequencies. It is clear that 
the precession frequency also increases with temperature 
as found for a particular value of r^, Fig. 0] (bottom). 

Up until this point in this paper we have focussed 
on systems with a fixed Nc, however, since experiments 
routinely conserve A^tot, all remaining simulations of 
this paper will be performed with a fixed total atom 
number, and variable condensate atom number. 

In the experiment of Freilich et al. (28j and in sub- 
sequent experimental runs by the group, the observed 
precession frequencies in a data set typically averaged 
a few (< 5) percent higher [7l| than those predicted 
for a zero temperature, Thomas-Fermi condensate in 
an axisymmetric trap 0, 0, [zS, IZ^]- next 
section we will continue our analysis on vortex precession 
frequencies for parameters of this experiment at both zero 
and finite temperature, in order to investigate the origin 
of this deviation. We will begin by briefly remarking on 
the procedure of the experiment before presenting our 
results for those parameters. 



A. Experimental Parameters of Freilich et. al [28| 

Vortices arise in this experiment during evaporation 
via the Kibble-Zurek mechanism [73 - [77j : the procedure 
for imaging is as follows: approximately 5% of the 
condensate atoms are outcoupled along the z-axis so 
that they are no longer confined by the trap, and they 
therefore, fall with gravity (along z). This proportion 
expands and the position of the vortex can be resolved 
and measured. This leaves the remaining 95% of the 
atoms trapped in the condensate and the vortex contin- 
ues to process in this slightly depleted condensate. At a 
later time, this process is repeated and the position of 
the vortex at that time also measured. The result is a 
series of images of vortex position which means that the 
real time dynamics of the vortex can be assessed. For 
the data presented in Ref. [l^ there was no discernible 
thermal cloud and the temperature was estimated to be 
T/Tc < 0.4. 

As a result of this technique, a single vortex was 
observed for approximately 655ms [78| . in a series of 
snapshots. Using these images, the precession frequency 
of the vortex line could be measured and the observed 
frequencies were found to average, in a typical data set, 
a few percent (< 5%) higher than those expected for 
the geometry and condensate parameters, given by the 
following equation; 

which uses the Thomas-Fermi approximation for the 
shape of the condensate. This experiment was conducted 
in a disk shaped condensate with a total atom number of 
Ntot « 4 - 6 X 10^ *^'^Rb atoms (we use Ntot = 6 x 10^ 
for our simulations). The trapping frequencies are as 
follows: ijj± = 27r X 36 Hz, and oj^ = Awj^, where A ~ a/S, 
therefore, the aspect ratio is similar to that used in the 
previous section, ensuring that the vortex stays relatively 
straight throughout its motion. 

Fig. [6] (top) shows the ZNG results (dashed lines) for 
precession frequency as a function of vortex position for 
various temperatures (increasing from bottom to top) 
and the predictions of Eq. @ for each temperature (solid 
lines). The bottom part of this plot gives the relative 
difference between the ZNG results for these tempera- 
tures and the corresponding zero temperature prediction 
obtained by solving the Gross-Pitaevskii equation. We 
solve the GPE for the number of atoms in condensate 
according to ZNG. We note that if we were to solve 
the GPE for the total number of atoms, the discrepancy 
between these results would be over 15% for T/Tc ^ 0.8. 

For a large condensate, it is natural to approxi- 
mate using Thomas-Fermi theory via the expression of 
Eq. ([3]), however, comparison of the zero temperature 
GPE results with the prediction of Eq. ^ , reveals that 
this prediction consistently underestimates the value of 
precession frequency. We find the relative difference 
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FIG. 6. (Color online) Top; Precession frequencies of 
a vortex for increasing temperature from bottom to top, 
T/Tc = (green diamonds, dashed line), 0.2 (39nK) (black 
triangles, dashed line), 0.4 (78nK) (red squares, dashed 
line), 0.7 (133nK) (blue circles, dashed line), and the 
respective analytically predicted frequencies from Eq. Q 
(correspondingly colored, solid lines). Bottom: Relative 
difference between the ZNG results and the GPE prediction 
at the vortex radial coordinate 0.25-Rtf. For completeness 
we indicate the values of Nc/N (top axis) in a nonlinear 
scale. The GPE results are obtained using the same number of 
condensate atoms as in the ZNG results. Trapping parameters 
as in text for iVxcT ~ 6 x 10^ *'^Rb atoms at Tc = 190nK with 
Rtf in the range « (9 — 9.5)/^ for the range of temperatures 
indicated. 



between these results to be in the range (5 — 10)% for 
the parameters shown. 

With increasing temperature, the ZNG results increas- 
ingly deviate from those obtained using the GPE for 
the same number of condensate atoms, therefore, the 
difference between the ZNG results and the predictions of 
Eq. ([3]) also increases with temperature (this behaviour 
is in roug h agreement with the experimental observa- 
tions [7l[), and is approximately 35% at the highest 
temperature of T/Tc — 0.7. Clearly, the TF formula 
(Eq. ([3])) is inherently approximate and should be used 
with care. 

A summary of the results of our simulations for the 
experimental parameters is shown in Fig. [T) The data 
for the decay rate, 7, and precession frequency, ujy, are 
all obtained with an initial vortex radial position tq = 
0.35i?TF- 
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FIG. 7. (Color online) Decay rate (left) and precession 
frequency (right) of a vortex at a position r„ = 0.35-Rtf, 
at temperatures T/Tc ~ 0.01 (green diamonds), 0.2 (39nK) 
(black triangles), 0.4 (78nK) (red squares), 0.6 (114nK) 
(magenta stars) and 0.7 (133nK) (blue circles). Parameters 
as in Fig. |6] 



IV. VORTEX CORE BRIGHTNESS 

The temperature of a system of bosons is commonly 
extracted by fitting a Gaussian profile to the high energy 
tails of the velocity distribution of the atoms [T^ • These 
atoms are associated with the high energy, noncondensed 
atoms in the system and can be described by a Maxwell- 
Boltzmann velocity distribution. A limitation of this 
procedure is that it becomes difficult to accurately 
extract the temperature at low temperatures when the 
thermal cloud density is low. Since the thermal cloud 
density is relatively high in the region of the vortex core, 
a vortex acts as a 'thermal-atom concentration pit' [sol ]. 
This led Coddington et al [s^l to suggest that measuring 
the density in the core region could be a possible tool for 
determining the temperature of the system. 

To perform a quantitative analysis of the thermal cloud 
density in the core, termed vortex core brightness [soj . 
we revert, for computational speed and efficiency, to a 
smaller system of atoms with a fixed total atom number 
of TVtot = 10, 000 ^■^Rb atoms. Fig. H] shows the 
condensate and thermal cloud local densities (left), and 
the densities integrated along the z— direction (right) in 
a harmonic trap with trapping frequencies uji_ = 2t: x 
129Hz and uj, — \/8lo±, at a temperature of T/Tc = 0.7. 
For both of these quantities, the vortex core appears as a 
dip in the condensate density with a corresponding peak 
in the thermal cloud density. An interesting property of 
the integrated thermal cloud profile is that throughout 
the extent of the condensate, the integrated thermal 
cloud density is relatively uniform, except for the peak 
that still emerges in the region of the vortex core. 

In Ref. [5^, the vortex core brightness was defined as 



B 



^cloud 



(4) 



where ricoro is the observed atom density, integrated 
along the line of sight at the core centre (see Fig. [5] 
right) and ridoud is the projected integrated density at 
this point based on a smooth fit of the overall atom 
cloud. We extract the integrated thermal density, ricore, 
at the central point of the vortex following this during its 
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FIG. 8. (Color online) Local density (left) and density 
integrated along the 2— direction (right), of the condensate 
(red) and thermal cloud (blue). The dip in the condensate 
density occurs at the position of the vortex. Trapping 
parameters as in text for A^tot ~ 10, 000 **^Rb atoms at a 
temperature of T/Tc = 0.7 where Tc = 177nK. 

motion, which we refer to as n|Q(rv, t). We calculate the 
projected integrated density, ndoud, at the same point as 
follows; we first run a simulation for the same parameters 
in the absence of a vortex and extract the integrated 
condensate density, henceforth denoted by ^^^(rv,^), at 
the given position. The value of ndoud is then obtained 
by the addition of the integrated thermal cloud density in 
the vortex core, n^]3(rv,t), with the projected integrated 
condensate density at that point, ri,^j-,(rv, t), I.e. ^cloud — 

''^2D + '^2D- 

We use the trapping parameters defined in Sec. IIIBI 
and above, noting here again that we fix the total number 
of atoms in the system to be consistent with routine 
experimental procedures. 

In Fig. [2] we plot vortex brightness, B, at the vortex 
position ry ~ 0.2i?TF for various temperatures. We see 
a clear trend in increasing brightness with increasing 
temperature for all approximations of the QBE, however, 
the value of the brightness is largest when the full QBE 
solved. We note again the importance of solving the 
full thermal cloud collisional dynamics in order to avoid 
underestimating the effect of the thermal cloud on vortex 
properties. 

More detailed information is provided in Fig. (TU] 
This figure shows the dependence of the integrated 
condensate n^Q(rv,<) (top row, right) and thermal cloud 
"2D (''v , t) (bottom) densities along the vortex trajectory 
(top row, left). The initial vortex position is 0.15i?TF 
and the temperature is T/Tc = 0.5. We also display 
densities for two other, higher temperatures. There 
is a clear trend of increasing integrated thermal cloud 
density with temperature. On the other hand, as 
one would expect, the projected condensate density 
decreases with increasing temperature. Furthermore, 
as the vortex moves radially outwards, the integrated 
condensate density decreases whereas, the integrated 
thermal density stays relatively constant (except near the 
edge of the condensate, not shown here). 

The results presented here indicate that core bright- 
ness, B, may indeed be a good candidate for extracting 
temperature of a condensate containing a vortex. How- 
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FIG. 9. (Color online) Vortex core brightness, B, for a 
vortex at position r-„ — 0.2_Rtf, as a function of temperature. 
Results for different levels of approximation are indicated by 
(i) solid, black dots: all collision (Ci2,C22) processes, (ii) 
open, magenta boxes: particle-transferring (C12) collisions 
only, (iii) blue crosses: thermal-thermal (C22) collisions only, 
(iv) open, green circles: no collisions, and (v) red stars: static 
thermal cloud approximation. Trapping parameters as in text 
for A^TOT = 10,000 ^■'Rb atoms, Tc = 177nK. Rtf varies 
between (4.2 — 4.9)/x for these results. 




FIG. 10. (Color online) Top: Vortex trajectory (left), 
Ti, — {xv,yv) for condensate containing a vortex with initial 
position ro = 0.15-Rtf at T/Tc = 0.5. Integrated projected 
condensate density, n2r,{^v,t) (right). Bottom: Integrated 
thermal cloud density, nJn(rv,t), at the centre of the vortex 
core, as a function of radial coordinate, for the temperatures 
T/Tc = 0.5 (red), 0.6 (magenta) and0.7Tc (blue). Parameters 
as in Fig. [9l 



ever, in order to make the method quantitative, one 
would have the take into account the dependence of B 
on the radial position of the vortex. Since remains 
practically constant throughout the motion of the vortex 
through the condensate (Fig. [10] (bottom)), our analysis 
suggests that this quantity itself may provide a better 
measure of the temperature of the system. 

We therefore extract the value of n^j^ for a range of 
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temperatures in Fig. [TTJ For completeness, we also con- 
sider the effect of the various collisional approximations. 
We see that the sensitivity to the collision processes 
included is much less than that of the vortex decay rate 
or precession frequency (see Fig.|4]). 



the temperature than total vortex brightness previously 
proposed by Coddington et al. Phys. Rev. A 70 063607 
(2004). 
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FIG. 11. (Color online) Integrated thermal cloud density, n^u, 
at the centre of the vortex core as a function of temperature. 
Results for different levels of approximation are indicated by 
(i) solid, black dots: all collision (Ci2,C22) processes, (ii) 
open, magenta boxes: particle-transferring (C12) collisions 
only, (iii) blue crosses: thermal-thermal (C22) collisions only, 
(iv) open, green circles: no collisions, and (v) red stars: static 
thermal cloud approximation. Parameters as in Fig. [5] 



V. CONCLUSION 

In summary, we have investigated observable proper- 
ties of vortices at finite temperature, including decay 
rate, precession frequency and vortex core brightness, us- 
ing a model which accounts for all collision processes be- 
tween the atoms. Particle-exchanging collisions between 
the condensate and thermal cloud provide the dominant 
contribution, while thermal-thermal collisions further 
affect these particularly at temperatures in the region 
T > O.bTc- While the decay rate of a vortex is dependent 
on initial vortex position, the precession frequency of a 
vortex is instead a function of the instantaneous vortex 
position; both quantities increase with increasing tem- 
perature. Furthermore, we have investigated the vortex 
precession frequency for the experimental parameters of 
Freilich et al. Science 329, II82 (2010), and found that 
even at low temperatures, there is a deviation from the 
predictions of a Thomas-Fermi analysis. This deviation 
increases with temperature and can be as large as 35% 
and should be detectable in future experiments. 

We also found that the integrated thermal cloud 
density in the region of the vortex core remains relatively 
constant as the vortex spirals towards the edge of 
the condensate. This observation suggests that the 
integrated thermal density is perhaps a better probe of 
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Appendix A: The mean-field potential and collision 
integrals 

Within the ZNG formalism, the thermal excitations 
are assumed to be semiclassical moving in a HE potential, 
i.e. an excitation with momentum p possesses energy 
ii = /2m -\- [/off (r, t) where the effective potential, C/eff 
is defined as 

l/eff(r,t) = Fext(r) -K2.gK(r,t) +n(r,0]. (Al) 

The quantities C22 and C12 appearing in Eq. ([2]) are 
collision integrals. C22 describes the redistribution of 
thermal atoms as a result of collisions between thermal 
atoms while C12 describes the change in the phase-space 
distribution function /(p,r, t) as a result of particle- 
exchanging thermal-condensate collisions. These are 
respectively defined as 



C22[f] 



dp2 



dp3 



dp4 



J {2T:hf J (27r?i)3 
X {2Tlhf5{p + P2 - P3 - P4) 

X 5{e + £2 - £3 - £4) 

X [(/ + l)(/2 + l)/3/4-//2(/3 



(27r?i)3 



and 



in 



dp2 



l)(/4 + l)], 

(A2) 



dp4 



h'' J {2TThf J {2TThY 

X {2^:Tlf6{m^rc + P2 - P3 - P4) 

X (5(£c + £2 - £3 - £4) 

X {2TThf[5{p^p2)-5{v-V7d 

X [(/2 + l)/3/4-/2(/3 + l)(/4- 



{2TTnf 



- 5{p - P4)] 
1)]. (A3) 



The delta functions in these expressions enforce conser- 
vation of energy and momentum. In particular those in 
the C12 term take into account that the condensate atoms 
have energy = mv'^/2+^c, and momentum mVc, where 
fic is the condensate chemical potential. 

If the condensed and noncondcnscd components are 
in diffusive equilibrium, C12 = 0. When they are out of 
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equilibrium, this term acts to transfer atoms between the 
condensate and thermal cloud. This exchange results in 



the source term of Eq. ((T]), 

(A4) 

The effects of C12 in the kinetic equation together with 
R in the generalised GPE ensure that the total particle 
number is conserved. 
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